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Abstract: We study the quark number susceptibility, an indicator of QCD phase transi- 
OO ! tion, in the hard wah and soft wall models of hQCD. We find that the susceptibilities in 

both models are the same, jumping up at the deconfinement phase transition temperature. 
We also find that the diffusion constant in the soft wall model is enhanced compared to 
^ ■ the one in the hard wall model. 
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1. Introduction 



There has been much interest in applying the idea of AdS/CFT|]l| in strong interaction. 
After initial set up for N=4 SYM theory, confining theories were treated with IR cut off 
at the AdS spaceQ and quark flavors |^] were introduced by adding extra probe branes. 
More phenomenological models were also suggested to construct a holographic model dual 
to QCD |§, m, 0]. The finite temperature version of such model were suggested in |8|, |9|, 
P^. For the purose of the Regge trajectory, quadratic dilaton background was introduced 



Ref. 1 11 1 whose role is to prevent string going into the deep inside the IR region of AdS 
space by a potential barrier and as a consequence the particle spectrum rise more slowly 
compared with hard-wall cutoff. Remarkably such a dilaton-induced potential gives exactly 



the linear trajectory of the meson spectrum. In |12], it was argued that such a dilatonic 
potential could be motivated by instanton effects. 

The quark number susceptibility Xgj which measures the response of QCD to a change 
in the chemical potential, was proposed as a probe of the QCD chiral/deconfinement phase 
transition at zero chemical potential p3|, 14 1. It is one of the thermodynamic observables 



that can reveal a character of chiral phase transition. The lattice QCD calculation |13| 
showed the enhancement of the susceptibility around Tc by a factor 4 or 5. Since then, 
various model studies [15, 16, and lattice simulations |18, 1£, 20, 21, |2^, ^] have been 
performed to calculate the susceptibility. 
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In this work, we calculate the quark number susceptibility and study QCD chiral/deconfinement 
phase transition in holographic QCD models ||6|, ^, 0. In the models adopted in the 
present work, we implicitly assume that chiral symmetry restoration and deconfinement 
take place at the same critical temperature Tc. 

In a gravity dual of QCD-like model, the confinement to de-confinement phase tran- 
sition is described by the Hawking-Page transition (HPT). At low temperature, thermal 
AdS dominates the partition function, while at high temperature, AdS-black hole geometry 
dominates. This was first discovered in the finite volume boundary case in |^], and more 
recently it is shown in [^] that the same phenomena happen also for infinite boundary 
volume, if there is a finite scale along the fifth direction. In our work, both hard wall Q 
and soft wall |11| models are considered, where we have a definite IR scale, so that we are 
dealing with theories with deconfinement phase transition. 

In the presence of the AdS black hole, the most physical boundary condition is the 
infalling one. In the hard wall model Q that the infalling boundary condition in the 
zero frequency and momentum limit can be understood as a conformally invariant Dirichlet 
condition. In the the soft wall model |11|, we calculate the quark number susceptibility 
and find that it is the same with the one obtained in the hard wall model, while diffusion 
constant of the soft wall model is enhanced compared to the hard wall model. 

The rest of the paper goes as follows. In section 2, we briefiy summarize the holographic 
QCD models adopted in the present work and discuss the chiral symmetry restoration in 
the models. In section 3, we calculate the quark number susceptibility in the AdS black 
hole background, adopting the infalling boundary condition, in the hard wall and soft 
wall models. Section 4 gives summary. In Appendix A, we re-evaluate the susceptibility 
with the Dirichlet boundary condition followed by the implication of the Hawking-Page 
transition |25| 



2. Holographic QCD and chiral symmetry at finite temperature 
2.1 Chiral symmetry in Hard wall model: 

The action of the holographic QCD model suggested in |6|, |^, is given by 



Si = y (i^xdzy/gCz 



4^5 



(2.1) 



where Dm^ = — iLm^ + i^Rhi and Lm = L%jT°-/2 with being the Pauh matrix. 
The scalar field is defined by <I> = Se^'^'''^" and < S >= ^v{z), where S is a real scalar 
and vr is a pseudoscalar. Under SU(2)y, S and vr transform as singlet and triplet. In this 
model, the 5D AdS space is compactified such that < z < Zm- ^ 

^ Zm is a infrared (IR) cutoff, which is fixed by the rho- meson mass at zero temperature: Xjzm — 
320 MeV, and the value of the 5D gauge couphng g\ is identified as g\ = VI-k^ jNc through matching with 
QCD 1,00. 
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As in , we work on the 5D AdS-Schwarzchild background, which is known to describe 
the physics of the finite temperature in dual 4D field theory side, 



dsl = ^ (^fiz)dt' - (dx'f - J^^dz^) ' /(^) = 1 - ^' (2-2) 

where i = 1,2,3. Here the temperature is defined by T = 1/{ttzt), and the Hawking-Page 
transition [25| occurs at Tc = 2^^^/{TrZm)- For the temperature lower than Tc, thermal 



AdS dominates, and it is hard to find temperature dependence in low temperature regime. 



which is actually consistent with the earlier work on large N gauge theory |24]. 
The equation of motion for v{z) in the black hole background is 



' zf + z^f 



t;(z)=0, (2.3) 



and the solution is given by S 



v{z) = M,z i, i, ^) + ^F,{\, ^, I ^). (2.4) 

Here Mq and Eg are identified with the current quark mass and the chiral condensate 
respectively. Note that at z = zt, both terms in v{z) diverges logarithmically. This 
requires us to set both of them to be zero, 

Mq = 0, = 0. (2.5) 

The latter condition means that the chiral symmetry is restored, and deconfinement and 
the chiral phase transition take place at the same temperature. It is interesting to observe 
that the mass term is also forbidden in this phase. This is consistent with the fact that the 
chiral symmetry forbids fermion mass term. ^ In reality the chiral symmetry is partially 
broken in low temperature and also partially restored in the high temperature due to the 
current quark mass. In this sense, the hard wall model respects the chiral symmetry more 
than the reality. ^ 



2.2 Chiral symmetry in Soft wall model 



In [11 1, dilaton background was introduced for the Regge behavior of the spectrum. 

Sii = / d^xdze-*£5 , (2.6) 



where $ = cz^ . Here the role of the hard- wall IR cutoff Zm is replaced by a dilaton- induced 



^Notice that we never discussed any fermions in this formalism, and hence we never defined any chiral 
symmetry in this theory explicitly. 

^In D3-D7 system, Tjg, Mq are not necessarily zero even for deconfined phase. 
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potential. The equation of motion for v{z) is given by 



zf z^f 



v{z) 



0. 



(2.7) 



At zero temperature [^], where / = 1, one of the two hnearly independent solutions of Eq. 
(|2.7| ) diverges as z ^ oo, and so we have to discard this solution. Then chiral condensate 
is simply proportional to Mq Now we consider a finite temperature case. Near zt, c 
dependent term in Eq. ( |2.7| ) is negligible and there is no difference between hard and soft 
wall model near the horizon. So we can draw the same conclusion of the complete chiral 
symmetry restoration. 



3. Quark number susceptibility 



In this section, we calculate the quark number susceptibility at high temperature in decon- 
finement phase. The relevant background is the AdS black hole, 



1 



(3.1) 



where \i={z/zt)'^, f(u) = and T = l/vrzy. 

The quark number susceptibility was proposed as a probe of the QCD chiral phase 
transition at zero chemical potential |13, ^ 



Xq 



dUg 



(3.2) 



The quark number susceptibility can be written in terms of the retarded Green's function. 
Here we follow the procedure given in |15]. First we write 



Xq{T,^i)=P J d-'xGoo{0,x), 
where the vector correlator at finite temperature is defined by 

-l/T 



(3.3) 



ah 



dT / d^xe-'P-"" I^J^{T,x)J^{fd,^) 



After the Fourier transformation, we obtain 



/oo 
—Goo{u},k) 



(3.4) 



(3.5) 



*In this work, we will not distinguish flavor singlet and non-singlet susceptibilities, xs and Xns respec- 
tively. Lattice QCD studies showed that xs — Xns [isl , [22| , which is because the mixing between 
isospin singlet and triplet vector mesons, p and uj, is tiny ]15[7^ 



- 4 - 



By the fluctuation-dissipation theorem, we have 



Goo{uj,k) 



1 - e-^/T 

where GQQ{ui,k) is retarded Greens function of Then we arrive at 



(3.6) 



hm 



du) 



27r 1 - e-""/^ 



(3.7) 



The ImG^(u;, k) is related to the real part of Green's function through the Kramers-Kronig 
dispersion relations, 

duj' ImGoo(<^, k) 



-( 



vr w — w 



One can easily see that Xq is written in the following form, 

X,(T,/.g) = -limReG^o(0,A:). 

What is the physics we want to see? The lattice QCD calculations 
the enhancement of Xq around T^. 



(3.8) 



(3.9) 



, 18, 19, showed 



^ XqjTc + e) 
^ - Xq{Tc - e) 



5. 



(3.10) 



The enhancement in may be understood roughly as follows |14]. At low temperature, 
in confined phase Xq will pick up the Boltzmann factor q~^n/t ^ where Mjv is a typical 
hadron mass scale ~ 1 GeV, while at high temperature the factor will be given in terms 
of a quark mass e"^^'/"^, and therefore there could be some enhancement. In Ref. it 
is shown that the enhancement in Xq may be due to the vanishing or a sudden decrease of 
the interactions between quarks in the vector channel. 

In the holographic QCD models, due to the HPT |25], the quark number susceptibility 
is described by the AdS black hole background at high temperature and by the thermal 
AdS at low temperature. We will patch them together. We calculate the susceptibility 
both in the hard wall model and in the soft wall model and describe how it changes under 
the phase transition. 



3.1 Hard vi^all model 

The temperature dependence of confinement phase can not be extracted from the AdS 
black hole due to the Hawking-Page transition. That is, thermal AdS replaces the AdS 
black hole at low temperature. Therefore we need to consider two phase separately. 



3.1.1 Confined phase 

We take the thermal AdS with hard wall at the dual gravity background. From 

the quadratic part of the 5D action Eq. ( |2.1| ), we obtain the equation of motion for the 
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time component vector meson 



z 



Vb(^,g1=0. (3.11) 



Note that in the above equation, the ^ Q hmit is already taken, following the definition 
given in Eq. (|3.9| ). With a finite go the equation will couple to other components. We 
solve Eq. ( |3.11 ) in a limit where g — > 0. The solution is Vb = ai — 02^^ . According to the 



AdS/CFT dictionary, Vq is dual to the number density j^{x) in boundary, hence we can 
identify oi as the chemical potential /i and the 02 as the charge density Q. Therefore 

Vq = lx-Q{z/z^m? . (3.12) 
Then the quark number susceptibility is given by 

1 9,Vb 



. (3.13) 

55 



For the black hole background, we have to request Vq = at the horizon for the regularity 
of the vector field, making // proportional to Q. However, for the thermal AdS there is 
no such requirement to be imposed. Therefore, even when a charge is zero, one can still 
have a finite chemical potential in confining case. However, in the limit where the chemical 
potential is zero, the charge density must be zero, otherwise we would need work to bring 
a particle into the system, requesting a finite chemical potential. That is, in the limit of 
zero chemical potential we should have zero charge and hence 

Xq = 0, (3.14) 

in the confined phase. 
3.1.2 Deconfined phase 

Here we calculate the quark number susceptibility in deconfined phase on the metric (p.lj). 
Note that the role of the IR cutoff z^ is none in the deconfined phase. The action of the 
gauge field, which is dual to the 4D quark current = q{t, x)^ij,q{t, x), is 

S = - [ d^x^^FMNF^'"" , (3.15) 
where I/55 = NcNj /{27r)'^ coming from the D3/D7 model |3(;]. In the presence of the 



AdS black hole, the most natural boundary condition is the infalling boundary condition 
since the black hole can only absorb classically. For the static case, we may consider the 
Dirichlet and Neumann boundary conditions. To impose the infalling boundary condition, 
we solve the problem at small but non-zero frequency and momentum and then take them 
to be zero. This is the problem considered in the literature on hydrodynamics [^]. Notice 
that in terms of the re-scaled coordinate u = (z/zt)"^, the momentum k = {to, 0, 0, q) enters 

2wT ' ^ 27rT 



in the action only in the combination of to = ly^, <\ = The relevant equations of 
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motion for the vector fields in the Au = gauge are 

f 



(3.16) 
(3.17) 



where ' means du- Out of the coupled equations we can eliminate A3 to get a second order 
differential equation for A'q := ^. 



{ufY . i»2-q2/(ti), , 



uf 



up 



(3.18) 



After taking out the near horizon behavior ~ (l — u) ^"^/^ of infalling wave, one can extract 
the small frequency behavior of the residual part using Mathematica. The result is 



2u^ 



^ = (i-.)--/2.5!4±^(i + !^i,. 

^ ^ (it»-q2) V 2 1 + n 



In 



2u 
1 + u 



+ 



(3.19) 



where A^ is the boundary value of Aii. With the prescription for the retarded Green 



function 



with 



6^S 



9i 



5Al{-k)5Al{ky 
^{-k)^{k) + • ■ 



one can get the retarded Green function 



G^,{u,q) = -i / dSe-^«-0(t)([J^(x), urn 



in small frequency; 



= 'y(.u,-,^Wln2/2) (1 + 1" 2 (r - 



(3.20) 
(3.21) 

(3.22) 
(3.23) 



Finally, by using eq. (|3.9D , we get the susceptibility 



(3.24) 



Xg/r^ is a constant which is consistent with high temperature behavior of lattice result 
P^ , [l8| , p^ , 2C, 21, 23]. Our results show that Xq{T)/T'^ jumps from zero, correspond- 
ing to the thermal AdS phase, to a constant, AdS black hole phase, with increasing the 
temperature. We plot Xq ^ ^ function of temperature. 



We close this section with an interesting observation. If we impose the Dirichlet con- 
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dition at the horizon, Vo{zt) = h, we get 



27r2 

X,(.T) = —{l-h)T\ (3.25) 

Notice that for the conformally invariant choice h = 0, Dirichlet condition gives the identical 
result to the infalling BC obtained above. 



3.2 Soft wall model 

Now we consider the soft wall model. The action is given by 

S = - [ d'x^e-'^^FMNF'''', $ = ^u. (3.26) 



tT) 

In this model, there is also the Hawking-Page transition ||2^. The equation of motion for 
Vq in the static-low momentum limit reads 

d,i-e~''^"d,Vo) = 0, (3.27) 
z 

whose solution is given by 

~ ci + C2Z^ ,z^0. (3.28) 

Then, following the same argument given in ( |3.1.1| ), we conclude that the quark number 
susceptibility is also zero in the soft wall model at low temperature. 

Now we consider high temperature phase. The relevant equations of motion are 

At: A'^-j^A[-^{q^At + cwA,) = (3.29) 
^o.: K + i0r2 + iU^K + ^i^-<^')Aa = O (3.30) 

■■ ^'z + {j0yl + 7@)^'. + llJW^^^^t - x^^A,) = (3.31) 
Au : toA't + qf{u)A', = (3.32) 



From the eq.( ^.29 ), we can express in terms of At 



Inserting it to eq.( 3.32| ), we obtain 



where ct = c/(7rT)^. Imposing the infalling boundary condition at the horizon, we take 
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4 = (1 - F^^j^ to obtain 



^// / c 1 — S-u^ w , 
[ttTY uj{u) I — u 



c 



l-3n2 /i(l + 2u) c i(l + n) 



+ V (vrT)2^ ^^"'V 2u/(n) (7rT)2 2f{u) 

W ^-//+;)^ -4-)f = 0. (3.35) 

In the long-wavelength, low frequency limit, F(u) can be expanded as series in t», q^, 

F{u) = Fo + voFi (u) + q^Gi + . . . (3.36) 
After some algebra, we obtain first a few terms 

Fo = B (3.37) 



Fi = 



12 



^2(e^^ - e"^") + (-12 + 5cT)e"^(^+") [Ei{-CT) - Ei{-CTu))) 
_2e'^T(2+w) ^ CT){Ei{-2cT) - Ei{-ct{1 + n))}^ (3.38) 
Gi = Se'^T" (^£^i(-cr) - Ei{-CTu) + e"^ (^Ei{-CT{1 + u)) - Ei{-2cT)) ^ , (3.39) 

where Ei(x) is defined as the principal value of Ei{z) = — j tdt. The integration 
constants of Fi, G\ are chosen by the regularity condition at u=l. B is determined from 
the boundary values of At and at u=0. 

B - f'rf (3.40, 
ixx){\ - ^cr) - 

Prom these results, J^tiu) is 

A!lu) = (1 - ur-l \ ^"fX^^f^ U + q^X,. + ^^Y^) , (3.41) 
ivo{l - ^ct) - q^ V 12 y 

where X^2 , Y^jj are 

' (^Ei{-CT) - Ei{-CTu) - e"^ (^Ei{-2cT) - Fi(-CT(1 + «))) j 
' ^2e"^^" - 2e-^^ + (-12 + 5ct) (fz(-ct) - Fi(-CT«)) 

-2e^^(-3 + ct) (Fi(-2cT) -Fi(-cr(l + 'ti)))^ • (3.42) 
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The real part of the retarded Green's function 

R<o(fc) = - g2 p2^2%2g4 [^W + D\^X^.{e) - -^D^PYroie)) , (3.43) 

where D = 1/2'kT and P = 1 — 5c/127r^r^. Note that the Green function has diffusion 
pole, and the diffusion constant is 

D = — \ . (3.44) 

i - 12;;^ 

Notice that it is dressed by the factor 1/P due to the effect of the soft wall. One can easily 
check the positivity of the diffusion constant in the relevant temperature regime 

P = 1 -0.17-1 >0, \iT>T^ (3.45) 

where Tc = l/vr^Zc^, and we used czc^ = 0.42 |2|]. 

The quark number susceptibility is obtained with eq.(|3.9|) and we get 



27r2T2 ^ ^ 

Xq = , 3.46 

which is the same with the result of the hard wall model. 

Putting together the results at low and high temperatures, we arrive at the following 
conclusion. Xq is zero up to the phase transition temperature, and it jumps to a finite 
value given in Eq. (3.46), which implies a first order phase transition between low- and 



high-temperature phases. The sharp transition might be the large Nc artifact. 

In hard wall case we observed that the infalling and (a specially chosen) Dirichlet 
boundary conditions give the same results. One may wonder if one can arrive at the same 
conclusion in the soft wall model. In Appendix we dig into this question to observe that 
those two boundary conditions lead to different susceptibilities. 

4. Summary 

We first discussed the chiral symmetry restoration in AdS/QCD models. The AdS/QCD 
models respect the chiral symmetry more rigidly than the reality in the sense that, in chiral 
symmetry restored phase, both of the chiral condensate and the mass of the quarks are 
zero. 

Then, we calculated the quark number susceptibility in both hard wall and soft wall 
models. At low temperature, in confined phase, we showed that Xgi which is defined 
in the limit of zero chemical potential, is zero. With the infalling boundary condition, we 
could uniquely determine the overall normalization of the susceptibility, unlike the Dirichlet 
boundary condition. We found that the susceptibilities in both models are the same with 
the infalling boundary condition, and ~ af high temperature, which is consistent 



with high-temperature lattice QCD observations [|^, 19, 21, 22, 23]. 
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In Appendix A, we considered Dirichlet boundary condition in the soft wall model. 
With the HPT, we predicted the temperature dependence of Xq &t high temperature apart 
from the overall normalization that is fixed by an IR boundary condition. Our result with 
the HPT exhibits a similar behavior observed in model studies [15, 16, 17| and lattice 
simulations |15, 19, 



21 



Regardless of the IR boundary conditions, our results in both models predicted a sharp 
jump in the quark number susceptibility, which is an unavoidable aspect of the HPT and 
could be smoothed out by including large Nc corrections. 

Finally, we discuss a limitation of our approach in the light of the QCD phase transition. 
The nature of the QCD transition depends on the number of quark flavors and the quark 
mass: for pure SU(3) gauge theory, it is a first order, for two massless quarks, it is a second 
order, for two quarks with finite masses, it is a cross over, for three degenerate massless 
quarks, it is a first order, etc. Unlike the Polyakov loop or chiral condensate, the quark 
number susceptibility is not an order parameter, and so in the present study we are not 
able to determine the order of the QCD phase transition. The susceptibility could serve, 
at best, as an indicator of the transition. 



A. Dirichlet boundary condition in soft wall model 

Here we give analysis with Dirichlet boundary conditions in high temperature. The equa- 
tion of motion for Vq in the static, zero momentum limit is the same as eq.( t3.27 ) and the 



2 

solution is still given by Vb = ae'^^ + h. However, with the Dirichlet boundary conditions 
Vo(0) = 1, Vq{zt) = 0, the result is given by 

We note here that Xq{To) ~ 1-22^^, where Tq = ^/c/^: is the temperature scale generated 
by c. 



B. parameters of the soft- wall models 

The masses of the vector mesons are given by = 4c(n + 1). If we use nii = 770 MeV 
and m2 = 1450 MeV to calculate the slope of the Regge trajectory, then we obtain ^/c ~ 
614 MeV and so end up with the reasonable value of the transition temperature Tc — 
195 MeV. The value of Tn was determined in Ref. |3l|, Tc = 210 MeV, and more recently 



the relation between Tc and Zm (\/c) is obtained through Hawking-Page analysis in the 



holographic models used in this work, where Tc ~ 191 MeV [25|. We note here that in 



1 32], the value of \/c was determined to be ~ 671 MeV. Finally, we relate c with the QCD 
string tension a. The masses of vector towers are given, in terms of cr, by = 2TTan. 
From this and = 4c(n + 1), we get c = ^o", so that the dilaton factor becomes e~^'^^ . 
Note that the relation between c and string tension was also observed in Ref. 
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(a) (b) 

Figure 1: (a) Xq in hQCD with tlie innfalling boundary condition discussed in section 3. (b) x? 
in the soft wall model with the Dirichlet boundary condition. The circles are for the (quenched) 
lattice QCD results as shown in Ref. |l^ , and the triangles are for full lattice QCD in |Q . 
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